Half-life of the Doubly-magic r-Process Nucleus 78 Ni 
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Nuclei with magic numbers serve as important benchmarks in nuclear theory. In addition, neutron- 
rich nuclei play an important role in the astrophysical rapid neutron-capture process (r-process). 78 Ni 
is the only doubly-magic nucleus that is also an important waiting point in the r-process, and serves 
as a major bottleneck in the synthesis of heavier elements. The half-life of 78 Ni has been experimen- 
tally deduced for the first time at the Coupled Cyclotron Facility of the National Superconducting 
Cyclotron Laboratory at Michigan State University, and was found to be HOlgg ms. In the same 
experiment, a first half-life was deduced for 77 Ni of 128^33 ms, and more precise half-lives were 
deduced for 75 Ni and 76 Ni of 344jJ° ms and 238±i| ms respectively. 



Doubly-magic nuclei with completely filled proton and 
neutron shells are of fundamental interest in nuclear 
physics. The simplified structure of these nuclei and their 
direct neighbors allows one to benchmark key ingredi- 
ents in nuclear structure theories such as single-particle 
energies and effective interactions. Doubly-magic nuclei 
also serve as cores for shell model calculations, dramati- 
cally truncating the model space, thus rendering feasible 
shell model calculations in heavy nuclei. All this is of 
particular importance for nuclei far from stability, where 
doubly-magic nuclei serve as beachheads in the unknown 
territory of the chart of nuclides 0, Q . 

When considering the classic nuclear shell gaps and 
excluding superheavy nuclei, there are only 10 doubly- 
magic nuclei, and only four of these are far from stabil- 
ity: 48 Ni, 78 Ni, 100 Sn, and 132 Sn. Of these, 48 Ni and 
78 Ni are the most exotic ones, and the last ones with 
experimentally unknown properties. 78 Ni therefore rep- 
resents a unique stepping stone towards the physics of ex- 
tremely neutron-rich nuclei. In a pioneering experiment, 
Engelmann et al. were able to identify three 78 Ni 
events produced by in-flight fission of a uranium beam at 
the Gesellschaft fur Schwerionenforschung (GSI), demon- 
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strating the existence of this nuclide. We report the first 
measurement of the half-life of 78 Ni at Michigan State 
University's National Superconducting Cyclotron Lab- 
oratory (NSCL), demonstrating that experiments with 
78 Ni are finally feasible. Such a measurement provides a 
first constraint for nuclear models and can serve as a first 
indicator of nuclear properties far from stability (See for 
example 0.). 

Very neutron-rich nuclei play an important role in the 
astrophysical rapid neutron-capture process (r-process) 
0, Q . The r-process is responsible for the origin of about 
half of the heavy elements beyond iron in nature, yet its 
site and exact mechanism are still unknown. 78 Ni is the 
only doubly-magic nucleus that represents an important 
waiting point in the path of the r-process, where the re- 
action sequence halts to wait for the decay of the nucleus 
0. 

One popular astrophysical site for the r-process is the 
neutrino driven wind off a hot, newborn neutron star in a 
core-collapse supernova explosion ||. In this case the r- 
process begins around mass number A = 90, with lighter 
nuclei being produced as less neutron-rich species in an 
a-rich freeze-out. For such a scenario Ni would not be 
directly relevant. However, the cv-rich freezeout fails to 
accurately reproduce the observed abundances for nuclei 
with A = 80—90 0, and the associated r-process does 
not produce sufficient amounts of the heaviest r-process 
nuclei around A =195 [Io| . 

78 Ni is among the important r-process waiting points 
in models that try to address these issues. Examples 
include models that assume nonstandard neutron star 
masses ^lj , or that are based on a supernova triggered by 
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the collapse of an ONeMg core in an intermediate mass 
star 0. In these models the neutron-capture process 
begins at lighter nuclei and the half-life of 78 Ni becomes 
a direct input. Together with the other already known 
waiting points, 79 Cu and 80 Zn, the half-life of 78 Ni sets 
the r-process timescale through the N = 50 bottleneck 
towards heavier elements, and also determines the for- 
mation and shape of the associated A = 80 abundance 
peak in the observed r-process element abundances. The 
A = 80 mass region has recently gained importance in 
light of new observations of the element abundances pro- 
duced by single (or very few) r-process events as pre- 
served in the spectra of old, very metal-poor stars in the 
Galactic halo. These observations point to the possibility 
of two different r-processes [lj, LUj being responsible for 
the origin of light r-process nuclei below A < 130. Only 
with accurate nuclear data, especially around 78 Ni - 80 Zn, 
will it be possible to disentangle the various contributions 
from neutron-capture processes in different astrophysical 
sites, and to interpret the data on neutron-capture ele- 
ments expected from the many new metal-poor stars to 
be identified in ongoing surveys [lij . 

In this experiment a secondary beam comprised of a 
mix of several neutron-rich nuclei around 78 Ni was pro- 
duced by fragmentation of a 140 MeV/nucleon 86 Kr + 
primary beam on a 376 mg/cm 2 Be target at the NSCL 
Coupled Cyclotron Facility. The average primary beam 
intensity was 15 pnA. Fragments were separated by the 
A1900 fragment separator |l|| operating with full mo- 
mentum acceptance. A position sensitive plastic scintil- 
lator at the dispersive intermediate focus was used to de- 
termine the momentum of each beam particle at typical 
rates of 10 5 /s. A 100.9 mg/cm 2 achromatic Al degrader 
was also placed at the intermediate focus of the separator 
to provide increased separation. 

Each nucleus in the secondary beam was individually 
identified in-flight by measuring energy loss and time of 
flight, together with the A1900 momentum measurement. 
The time of flight was measured between two scintillators 
separated by about 40 m: one located at the intermedi- 
ate image of the A1900 and the other located inside the 
experimental vault. The beam was stopped in a stack of 
Si detectors of the NSCL Beta Counting System (BCS) 
[r?j |. Energy loss was measured in the first two Si de- 
tectors, which were separated by a passive Al degrader 
of variable thickness. The degrader thickness was ad- 
justed to stop the nuclei in a 985 fim double-sided Si 
strip detector (DSSD). The DSSD was segmented into 
40 1 mm strips horizontally on one side, and vertically 
on the other, resulting in 1600 pixels. The beam was 
continuously implanted into the DSSD, which registered 
the time and position of each ion. The typical total im- 
plantation rate for the entire detector was under 0.1 per 
second. 

Using the dual-gain capability of the BCS electronics, 
the DSSD also registered the time and position of any 
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FIG. 1: Particle identification using energy loss vs. time of 
flight for a subset of the data. 



/3-decays following the implantation of a nucleus. This 
allowed the correlation of a decay event with a previously 
identified implanted nucleus. Additonal Si detectors in 
front and behind the DSSD were used to veto events from 
light particles in the secondary beam that can be simi- 
lar to /3-decay events. With this setup, the total /3-type 
event background rate associated with an implanted ion 
was typically less than 3 x 10~ 2 /s. Fig.Q]shows the par- 
ticle identification using energy loss vs. time of flight. A 
total of 11 78 Ni events were identified over a total beam 
time of 104 hours. Using a rough estimate for the in- 
tegrated beam current and a transmission for Ni frag- 
ments into the experimental vault of 65% calculated with 
the Monte-Carlo beam transport code MOCADI 0] , we 
obtain a rough estimate of the production cross section 
of 0.02±0.01 pb. This is lower than the estimated cross 
section from the EPAX formula of 4 pb 0] . 

Decay half-lives were determined using a maximum 
likelihood analysis, which has been used before in ex- 
periments with low statistics [20i . and in extreme cases 
with just 6 and 7 /3-decay events [il]. For this work, 
the formalism was modified to account for /3-delayed neu- 
tron emission. The method finds the decay constant that 
maximizes a likelihood function, which is the product of 
probability densities for three decay generations as well 
as background events, to produce the measured time se- 
quence of decay-type events following the implantation of 
a beam particle. The calculation requires knowledge of 
the /3-detection efficiency, background rate, daughter and 
granddaughter half-lives, including those reached by (3- 
delayed neutron emission, and branchings for /3-delayed 
neutron emission [P n ) for all relevant nuclei in the decay 
chain. 

For 75 Ni, 76 Ni, 77 Cu and 78 Cu the statistics were suf- 
ficient to determine the /3-detection efficiency by com- 
paring fitted decay curves with the total number of im- 




planted species of that isotope. The resulting efficiencies 
agree very well and range from 40% to 43% with no sys- 
tematic trends in the deviation. For 77 Ni and 78 Ni, an 
average efficiency of (42 ±1)% was adopted. The back- 
ground was determined for each run (typical duration 
of lh) and in each detector pixel by counting all decay 
events that occur outside of a 100 s time window fol- 
lowing an implantation. Because of the low implantation 
rate the background is constant over the 5 s time window 
used to correlate decays to an implantation. Experimen- 
tal P n values as well as daughter and grand-daughter 
half-lives used for the analysis were taken from [22j and 
1 2 .'"'J] when available. The experimentally unknown P n val- 
ues for the Ni isotopes were taken from detailed spheri- 
cal quasi-particle random-phase (QRPA) calculations for 
pure Gamow- Teller (GT) and GT with first-forbidden 
decay [24| and a number of different choices of single- 
particle potentials and mass model predictions. From 
this study, we derive an average uncertainty for the cal- 
culated P n values of about a factor of two. 

The statistical error of the derived decay half-lives is 
obtained directly from the maximum likelihood analysis. 
As sources of systematic errors we considered uncertain- 
ties in the P n values and daughter or grand-daughter half- 
lives, as well as uncertainties in background rate and de- 
tection efficiencies. The systematic uncertainties for the 
half-lives of 78 Ni, 77 Ni, 76 Ni, 75 Ni are, (in ms) ±f Q , ±£ x , 
-5, -e, respectively. The main contribution to the sys- 
tematic errors are uncertainties in the detector efficiency, 
and uncertainties in the parent P n values. In the case of 
78 Ni, we also took into account the possibility that one 
of the events is misidcntified. Given the very low number 
of events beyond 78 Ni in the particle identification (see 
Fig. ^) this is a very conservative assumption. For 78 Ni 
this leads to a systematic error of lg°ms. 

Systematic and statistical errors are correlated since 
the shape of the likelihood function depends on the anal- 
ysis parameters. To add systematic and statistical errors 
we therefore reran the analysis for all combinations of 
systematic variations and employed the lower and upper 
one-sigma limits of the resulting statistical errors as the 
total error budget. 

In principle our analysis depends somewhat on the un- 
known feeding and decay branchings of the known iso- 
meric states in 76 Cu and 77 Zn, which are part of the 
decay chains considered here. Assuming decay from the 
isomeric state with a half-life of 1.27 s for 76 Cu would 
increase the 76 Ni half-life by no more than 12 ms and 
the 77 Ni half-life by no more than 5 ms. Assuming pop- 
ulation of the 1.05 s isomer for 77 Zn could change the 
half-life of 77 Ni by -8 ms to +13 ms, and the half-life of 
78 Ni by -10 ms to +15 ms depending on the probability 
for that state to /3-decay. These uncertainties are based 
on extreme assumptions with no obvious central value. 
We therefore give them separately and do not include 
them in our systematic error bars. 
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FIG. 2: Experimental Ni half-lives from this and previous 
work p2| compared to different theoretical calculations. 



For the assumption that Cu and Zn /3-decay from 
the ground states with half-lives of 0.641 s and 2.08 s 
respectively our final results are 3441^4 m s for 75 Ni, 
23S±^| ms for 76 Ni, 128t^ ms for 77 Ni, and 110±^° ms 
for 78 Ni. For 77 Cu and 78 Cu, we obtain 450l2i ms and 
3231 }g ms, in excellent agreement with previous work 
(469+8 ms and 342+11 ms) 25]. 

In Fig.|21our new experimental half- lives are compared 
with various theoretical predictions. Often employed in 
r-process model calculations are the global QRPA cal- 
culations of Moller et al. 1997 j2(| or Borzov et al. 
1997 H3, the latter being limited to spherical nuclei. 
Our results show that the trend of these models to over- 
predict Ni half-lives by factors of 3-4 already observed 
for the more stable isotopes persists into the path of 
the r-process at 78 Ni. The recent versions of both mod- 
els [23, besides other improvements now also include 
first-forbidden transitions. They clearly lead to better 
though still somewhat large half-life predictions. Fig. [21 
also shows results from calculations with the same model 
as Moller et al. 2003 [2^] but using a mass model that 
includes a quenching of shell gaps far from stability [ex- 
tended Thomas-Fermi approach + Strutinsky Integral, 
with shell quenching (ETFSI-Q) 30]]. These calculations 
give the best agreement with experimental data among 
the global models. 

In order to better understand the nuclear structure 
in this mass region and to benchmark global models 
beyond the range of experimental data it is important 
to test the more sophisticated microscopic calculations, 
which have been performed for a limited set of singly- 
and doubly-magic heavy nuclei (see Fig. |3J|. The self- 
consistent QRPA approach [3l| agrees with the shell- 
model calculation [32j and the experimental data for most 
nuclei, but predicts a 78 Ni half-life that even exceeds the 
experimental 80 Zn half-life. Our measurement clearly fa- 
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FIG. 3: Theoretical half-life calculations for N = 50 compared 
to experimental data from this and previous work [2^1 , Engel 
et al. only gives values for even Z nuclei. 



vors a much lower 78 Ni half-life. On the other hand, the 
shell model results are in good agreement with experi- 
mental data. Of course this does not necessarily mean 
that the shell-model description of this mass region is 
entirely correct. For example, deviations in excitation 
energies, transition strengths, and decay Q- value can in 
principle compensate each other. More experiments in- 
cluding detailed spectroscopy as they might become pos- 
sible at future facilities will be needed to clarify this. 

In summary, we present the first results for the half- 
life of 78 Ni and other neutron-rich Ni isotopes. With 
these results, experimental half-lives are available for all 
but one ( 48 Ni) classical doubly-magic nuclei. Also, the 
half-lives of all important N = 50 waiting points in the 
r-process are now known experimentally. This will make 
r-process model predictions of the nucleosynthesis around 
A = 80 more reliable and comparison with observational 
data more meaningful. It will also put the overall delay 
that the N = 50 mass region imposes on the r-process 
flow towards heavier elements on a more solid experi- 
mental basis. In this respect the half-life of 78 Ni is of 
special importance as during the initial stages of the r- 
process when the heavier nuclei are synthesized the r- 
process path passes through 78 Ni and 79 Cu rather than 
through the more stable N = 50 nuclei 33]. The de- 
lay timescale for the buildup of heavy elements beyond 
N = 50 is therefore set by the sum of the lifetimes of 
78 Ni and 79 Cu. Our experimental data clearly favor the 
short timescale of 450 ms obtained with the prediction of 
Langanke and Martinez-Pinedo |32| over the much longer 
delays of 960 ms predicted for example by Moller et al. 
|2rij | leading to an acceleration of the r-process. This is in 
line with recent improvements in theoretical /3-decay half- 
life predictions along the entire r-process path that also 
tend to result in shorter half-lives thereby speeding up 



the r-process |2jj. Detailed r-process model calculations 
with the new experimental data are beyond the scope of 
this paper, but will be presented in a forthcoming study. 
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